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Abstract 
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dent and N >2 dependent variables is considered. We choose the existence 
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Introduction 



In the papers [H El E] a general theory of integrable systems of PDEs of the form 

ut = A{u)u^ + B{u)uy, (0.1) 

where u is an A^-component column vector, A{u) and B{u) are A^x A^-matrices, was developed. 
The existence of sufficiently many of the hydrodynamic reductions |H H] has been proposed as 
the definition of integrability. Unfortunately, for arbitrary it is difficult to write down 
explicitly the conditions for A and B, which follows from this definition. Nevertheless, for 
N = 2 the complete set of integrability conditions have been found in the paper [2] . For N > 2 
even the verification whether a given equation is integrable is a serious task. Any attempt for 
classification of integrable models based directly on this definition seems to be hopeless. 

To overcome the difficulties for A^ > 2 the following two observations [21 [3] are very useful. 
First, under some conditions of generic position, the matrix M = (1 + tA)~^{l + kB) for 
integrable models should be diagonalizable by a point transformation u $(u) for generic 
values of the parameters t and k. If the eigen- values of M are distinct, this is equivalent to the 
fact that the Haantjes tensor [12] of M is identically zero. This gives rise to an overdetermined 
system of the first order PDEs for entries of A and B. Given a system (10. II) it is not difficult to 
verify whether these PDEs are satisfied or not. The simplest equations from this overdetermined 
system also can be useful for classification of integrable models (10.11) . 

The second observation made in [1] is that for N = 2 the integrability conditions are 
equivalent to the existence of the scalar pseudopotential 

^i = /($„u), =^(^„u), (0.2) 

for (larfl. The scalar pseudopotential plays an important role in the theory of the universal 
Whitham hierarchy [3 [6l [7j. A possible importance of pseudopotentials was also noticed in 
[8]. The existence of pseudopotential implies a representation of (10. ip as the commutativity 
condition for the corresponding characteristic vector fields 

d d d d „ d „ d . >.t 

For recent attempts to use similar representations for integration of dispersionless PDEs see 

mm- 

In this paper we assume that the matrix A{u) is constant and consider integrable systems 
of the form 

Uit = XiUi^+ ^ bij{u)ujy, (0.3) 



^The latter means that that the overdetermined system (|0.2p for is compatible if and only if u is a solution 
of (lO). 
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where Ai, . . . , Aat are pairwise distinct constants. Functions bij (as well as all other functions) 
are supposed to be locally analytic. Note that the transformation 

Ui Mui) (0.4) 

preserves the form of the system (10.31) for arbitrary functions of one variable i/jilui). 

For N = 2 such systems were considered in [2J. Our goal is to obtain a list of the most 
interesting examples of integrable models (10.31) with N > 2. As far as we know, nobody 
systematically investigated such systems before us. 

Example 1. Consider equation (10. 3p with 

_ Xj — \j . , . 

Ui — Uj 

bii = -^bij, (0.5) 

where cj are arbitrary constants. It is not difficult to verify that for arbitrary this equation 
possesses pseudopotential (10. 2p . where 

Af N 

g = '^Ci log{ui -"^y), f = Y1 ^°S(^ii - ^s/) • (0.6) 

i=l i=l 

This pseudopotential has the following structure: 

N N 
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J2hi{^,u,), / = ^A,/i,(e,n,), (0.7) 



i=l 



where ^ = "^y. In Section 2 we show that (10. 7p is true for any integrable equation of the form 

(ED. 

Notice that functions hi{^,Ui) in (10.61) have "moveable" singularities with respect to the 
variable ^ = "^y. This means that the position of the singularity depends on u. 

In this paper we describe all equations (10. 3p possessing pseudopotentials such that for any i 
the function hi{^, Ui) has a movable singularity. This leads to series of new interesting examples 
of integrable systems (10. 3p . These examples are presented in Section 1. It would be interesting 
to find the hydrodynamic reductions for these equations and describe the multiple waves [131 [I] 
in terms of these reductions. 

In Sections 2-4 we deduce a functional equation describing the pseudopotentials with move- 
able singularities and find all it's solutions. These solutions correspond to examples of Section 
1 and their degenerations. 
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1 Examples. 

Example 2. Consider equation (I0.3P given by 
and 

bit ^ ^ ^jj • 

Here if «: = — 1 or /c = 0, then Cj are arbitrary constants. For other k the constants Cj should 
satisfy the following two relations 

N N 

Y,Ci = 0, ^\c, = 0. (1.8) 

i=l i=l 



It is easy to verify that for any this equation admits a pseudopotential (10.71) with 

h,iC, u) = Q («:(e -u) + log(e«-" - 1)) . (1.9) 

Example 3. Consider equation (10. 3p with 

bij = Cj{Xj - Xi) _ ■ — , It J, 



where cj are arbitrary constants satisfying conditions (II. 8p . This equation has a pseudopotential 
(1071) with 

/i,(e, u) = c, ((e + 1) log(M) - log(« - O) • (1.10) 



The equation from Example 1 is a particular case of the following model: 
Example 4. Let 

_ (A, - Xj)cjP{ui)M{u,) . . 



Ui — Uj 

_ ^ ^ (A. - A,)c,P(n,)Af (..,) _ ^ _ _ 

^ Ui- Uj ^ 

where Cj are arbitrary constants, and the functions B, M are defined by quadratures from 
B' = ihx + h + z,) M, M' = MZh^±^l^±^. 
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Here 

P{x) = k^x^ + k2x'^ + kix + fco 
is an arbitrary polynomial of degree not greater than 3, zi, zq are arbitrary constants. 
The corresponding equation (I0.3P possesses pseudopotential (I0.7p . where 

, N 1 , ^ X P{u)M{u) 



u-m Mm))' ' ' ' u-m ' 

0' = P(0)M(0). 

For any given P and z\, zq the equations for (f){^) and it) can be easily solved by quadratures. 
Using admissible transformations 

auj + b . 
^ —7, t = l,...,N, 

cui + d 

one can reduce the polynomial P to a canonical form. For example, if all three roots of P are 
distinct, then without loss of generality we may put P{x) = x{x — 1). In this case 

M{x) =x-''{x-l)'\ 

where Si = —2:07 ^2 = Zq + Zi, and 

+ + 

It is not difficult to find that 



^(e, n) = ,,,,,,,,,,, / : dt - c^^+\c - ly 



0(e^(</>(o-i^7e m-t ' ' J m-c 

where 0' = (j)'^+\(f) - iy^+\ 

Other two canonical forms are P = x and P = 1. The latter generates Example 1 if 
zi = zo = 0. 

2 Pseudopotentials. 

A pair of equations of the form 

^t = fi'^y,Ui,...,UN), = gi^y,Ui,...,UN), (2.11) 

with respect to unknown \I' is called a pseudopotential for equation (10. 3p if the compatibility 
condition "^t^ = ^^.^ for (12.111) is equivalent to (10.31) . Differentiating (12. lip , we find that this 
compatibility condition is given by 

AT AT AT AT 

h ^i'^^)y 9i9 + ^{ut)x dif = ^iui)y dif + ^{ui)t d^g. 
1=1 1=1 1=1 1=1 
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Here and below we denote by ^ and ^ by di. Substituting the right hand side of flO.Sp 
for t-derivatives and sphtting with respect to x and ^/-derivatives, we get that for any i the 
following relations hold: 

dif = XAg, (2.12) 

N 

dig -g^dif = bjidjg. (2. 13) 

i=i 

Since Aj are pairwise distinct, it follows from the condition (12.121) that 

N 
i=l 

and 

N 
i=l 

It is easy to see that the integration constant c{^) can be distributed between functions hi. 
Thus we have arrived at (10.71) . 

Substituting (10. 7p into (I2.13p . we obtain 

dihi{^, Ui) ^{Xj - Xi)hj^{^, Uj) = ^ bji djhj{^, uj). (2.14) 
J j 

Remark 1. If we fix a generic value .^o of the variable ^, then it follows from (I2.14p that 

bu = Y.^\ - m^j) - E^^-*^' (2.15) 
where (t)j{uj) = hj^{^o,Uj), Sj{uj) = djhj{^o,Uj). 



3 Basic functional equation. 

Suppose that for any i the function /ij(^, u) has a singularity with respect to ^ and this singular- 
ity depends on u. After a transformation of the form (10.41) . we may assume that for each i the 
singularity of hi{^,u) is located on the diagonal ^ = u. We say that hi{^,u) has a singularity 
on the diagonal, if for fixed generic u and any e > we have max{|/ij((^, u)\,\^ — u\ < e} = oo. 

Proposition 1. Suppose hi{C,,u) has a singularity on the diagonal^ = u for each i. Then 
there exist: a function h{^,u), functions /j(0 ^^'^ non-zero constants Ci such that 

h,{C,u) = cMC,u) + fM), (3.16) 
bji = (Xi - Xj)ci dih{uj, Ui), i^ j, (3.17) 
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bii = "^i^j - K)icj dih{ui, Uj) + fjiu.i)), (3.18) 

h{^,u) = ln(^ — m) + regular part. (3.19) 

Moreover, the following functional equation 

h^{C, w) v) + K{w, v) hyj{i, w) - K{v, w) v) = z/(^, v) (3.20) 

holds for some function u. 

Proof. Considering (12.141) near the diagonal ^ = Uj, where j ^ i, and comparing the 
singularities, we obtain 

bji = (Xj - Xi)fij{uj) dihi{uj, Ui), (3.21) 

where fJ.j{uj) = lim^^„^. ^^0^-^- We see that the function hji depends on Wj, Uj only and has the 
same singularity on the diagonal as dihi{uj,Ui). Considering (I2.14p near the diagonal Ui = uj, 
comparing the singularities and using (I2.15p . we obtain djhj{^,Ui)si{ui) = dihi{^,Ui)sj{ui) or 
dihi{^,Ui) = h'i{ui)r{^,Ui) for some functions Ui and r. On the other hand, consider (12.140 for 
generic values ^i, ^at of variable ^. For each fixed z = 1, we have a system of linear 
equations for bn, ...,biiy with matrix Q = (qjk), where qjk = Vj{uj)r{^k,Uj). This system must 
have a unique solution by definition of pseudopotential and therefore A = det Q 7^ 0. It is clear 

p. . 

that bji = where Pji is regular on each diagonal Uk = ui. It is easy to prove the following 

Lemma. Let A(tii, m^) be the determinant of an m x m matrix Q, whose entries qij 
have the form qij = gi{uj) for some functions gi, ...,gm- The function A is not equal to zero 
identically iff the functions gi,...,gm are linearly independent. In this case diA 7^ on the 
diagonal Ui = uj for each i ^ j. 

From this lemma it follows that the only singularity of bij on the diagonal can be a pole of 
order one. Taking into account ( 13.210 . we obtain that near uj = Ui 

dihi{uj,Ui) = — ^-LjI. _)_ regular part 

Uj — Ui 

or, after integration, 

hi{uj,Ui) = —ai{uj) \n.{uj — Ui) + regular part. 

Considering the singular part of (12.140 at ^ = Uj, we obtain fJ-jiuj) = —1 and a[{uj) = 0, i.e. 
C(i{uj) = —Ci for some constant q. Comparing the singularities in (I2.14p at ^ = -Uj, we find that 

^ii = ^i^j - ^i) dihj{ui, Uj). 

Substituting this expression for ba into (I2.14p . we obtain 
^(Aj-Ai) (jij^{^,Uj) dihi{^,Ui) + dihi{uj,Ui) djhj{^,Uj)-dihj{ui,Uj) dihi{^,Ui)^ =0. (3.22) 
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Considering the singular part of (13.221) at Ui = Uj, we get 

Cj{hi{^,u))u = Ci{hj{^,u))u, 

which gives f l3.16p and f l3.19p . Now (13.171) and (13.181) follow from (I3.16P and expressions for bij 
and bii already obtained. Using (13.221) . we arrive at the relation 

Y.j=iii^j - {h^{^, Uj) dih{^, Ui) + dih{uj, Ui)djh{^, Uj) - dih{ui, Uj)dih{^, Ui) + 

(/;(0-/j(«.))We,n,)) =0. 

This relation gives the equation (I3.20p for some function i/{^,v). ■ 
Remark 2. If a pair h{^,v), i^{^,v) is a solution of (13.201) . then 

He, v) = h{^, v) + /(O, v) = v) + ifiO - f'{v))K{C, v) (3.23) 

is also a solution of (I3.20p . Therefore, if ^{C,, v) has a form {g{0^9i''^))^v{C,, v) for some function 
g, then we can bring v) to adding to h{^, v) a suitable function of ^. 

Remark 3. Without loss of generality, we can assume that fi{^) = fj{C,) for all Indeed, 
only the linear combinations X]i/*(0 ^^"^ appear in (10. 7p . Furthermore, according 

to Remark 2, we may put /j(0 = 0. 

Remark 4. It follows from (I3.17p . (I3.20p that for any i j the function b = bij{ui,Uj) 
satisfies the following functional equation 

b{w, v)buj{x, w) — b{x, v)by{v, w) + b{x, w)buj{w, v) + b{x, v)bx{x, w) = 0. 

Remark 5. It follows from (I3.17p . (I3.18P that the equation (10. 3p defined by (I3.17P and 
(I3.18P with fi{^) = can be written in the following divergent form: 

where 

(Ti = ^(Aj - Xi)cjh{ui,Uj). 

Thus any such equation has at least linearly independent hydrodynamic conservation laws. 

Proposition 2. Let h{C,,v) be a solution of (13.200 with z/(^,f) = 0, then for any nonzero 
constants Ci the formula 

N N 

g = ^Cih{^,Ui), / = X] ^i^Mi^ Ui) (3.24) 

1=1 i=l 

defines a pseudopotential for equation (10. 3p given by (13.170 and (I3.18P with fi{^) = 0. 
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Let h{^,v) be a solution of (13.201) with 7^ 0, then fl3.24p defines pseudopotential for 

equation (10. 3p given by fl3.17p and fl3.18p with fi{^) = iff the constants Ci satisfy the relations 

(O). 

Proof. According to the previous results, (I3.24p defines a pseudopotential iff 
y^(A-,- - \i)cj (ji^iC, Uj) dih{^, Ui) + dih{uj, Ui) djh{^, uj) - dih{ui, uj) dih{^, Uij^ = 0. 

Substituting (13.20p into this relation, we obtain the statement of the proposition. ■ 



4 Classification of solutions for the functional equation. 

Proposition 3. Let a pairh{^,v), i'{C,,v) be a solution of fl3.20p with asymptotic (13.190 . Then 
up to substitutions of the form fl3.23p it belongs to the following list: 

h{x, v) = k{x — v) + \og{x — v), iy{x, v) = k{k + 1); 

h{x, v) = k{x — v) + log(e'^~^ — 1), v) = k{k + 1), 

where n is an arbitrary constant] 

X 

h(x, v) = (x + 1) log(f ) — \og(v — x), fix, v) = 

V 

, r PU{x))(t)'{tf , [ (j)'(c) , , , 

h(x,v] = / , , , , , / XX -r^/ , , XX ,,/ . at ~ / — — -r-rdx, i'(x,v] = 0. 

Here c is a constant and the function (p is defined by the following differential equation: 

^ V3P((/)) P{(P)J ^ ' 

where 

P{x) = k^x^ + k2x'^ + kix + /co, Z{x) = zix + zq 
are arbitrary polynomials such that deg P < 3 and deg Z < 1. 
Proof. According to f l3.19p . we have an expansion of the form 

00 

h{w, v) = \n{w + ai{w){w - v)' (4.25) 

as V tends to w. To describe the solutions of the functional equation ( 13.20p . let us investigate a 
set of conditions relating the functions a^. Using expansion f l4.25p for h{v,w) and h{w,v) and 
equating the coefficients of different powers oi v — w, we obtain an infinite sequence of PDEs 
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for the function h{x,v) and the coefficients ai{v). The simplest three of these PDEs read as 
follows: 

- 2hyhr,^ + 2aiK^ = 0, (4.26) 
hvvvv — ^Khxvv + 3ai/it,t,t, + 6(% + 2a2)hyy + 3(a'/ + + I2a3)hy = 0, (4.27) 
hvvvvv ~ 4:h^hxvvv + ^ittihyyyy + 12(a'^ + 2a2)h^vv + 12(a'/ + Aa'2 + 603)/!^^ 

(4.28) 

+4(a;" + 64' + 124)/i^ = 0. 

Substituting expansion fl4.25p for /i(x, v) to (14.261) . we observe that all coefficients a^, i > 2 can 
be expressed as certain differential polynomials of ai and 02. For example, 

as = + 2aia'i + 4^). 

This means that the function h{x,v) is uniquely determined by functions ai{v) and a2(^)- 

The expansion of equations (14.271) and (I4.28P leads to differential relations between ai and 
02. In particular, the simplest relation following from (I4.27P has the form 

a'l' + Gaia'l + 6a[ ^ - 6ala[ + 12{a2a[ + a'^ai) = 0. 

If tti 7^ 0, this implies 

C - < - 6aia[ + 2al 
12^^ ^^-^^^ 

for some constant C. Eliminating 02 with the help of (I4.29p . we arrive at an overdetermined 
system of ODEs for function Oi. If C = 0, we find from this system that 

aj aS^^ + 2ai{3al - a[) af ^ - 4ai (a'/)^ + 2{a[ ^ - 9ala[ + Aa^) a'/ - 16a?a; ^ = 0. (4.30) 

The general solution of this equation can be written as 

3Z(^ f2n^ Zi^)^ 

"^^^^ = -2P(0M)'^(^^' = I 3P(0) + P(0)r ' 

where 

P{x) = k^x^ + k2x'^ + kix + fco) '^(^) = ^1^ + ^0 

are arbitrary polynomials such that deg P < 3 and deg Z < 1. For any given P and Z 
the equation for <p can be easily integrated by quadratures. Example 4 from Section 1 after 
transformation u — > 0(tt) describes the pseudopotential generated by such a function ai. 

Let C 7^ 0. Then a simple analysis of the ODE system for ai shows that either a'^ = or 

4:{a\f + 12a? {a[f + I2{a{ - Ca^) a\ + 4a? + 4Ca? + = 0. (4.31) 

It is easy to verify that if a'^ = 0, then a^ = for any i and therefore h{x, v) = H{x — v) for 
some function H. Solving equation (I4.26p . we get 

H{x) = ci + C2X + log(l - e'"^'') 
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or 



H{x) = Ci + C2X + log(c32;). 



These solutions correspond to the model of Example 2 and it's degeneration. 

The left hand side of fl4.3ip can be decomposed into three factors. Each factor gives rise to 
a differential equation of the form 



where k is related to the constant C from (14.291) by C = 2k^. The corresponding model is 
described in Example 3. 

The case ai = should be considered separately. It is easy to get that in this case 



It turns out that it is a particular case of the model described by f l4.30p . In the corresponding 
formulas one has to put Z = 0. M 
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